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Abstract. Effective Field Theory (EFT) provides a powerful framework that exploits a separation of scales in physical 
systems to perform systematically improvable, model-independent calculations. Particularly interesting are few-body systems 
with short-range interactions and large two-body scattering length. Such systems display remarkable universal features. In 
systems with more than two particles, a three-body force with limit cycle behavior is required for consistent renormalization 
already at leading order. We will review this EFT and some of its applications in the physics of cold atoms and nuclear physics. 
In particular, we will discuss the possibility of an infrared limit cycle in QCD. Recent extensions of the EFT approach to the 
four-body system and A'-boson droplets in two spatial dimensions will also be addressed. 
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INTRODUCTION 

A separation of scales in a physical system can be exploited using the framework of Effective Field Theory (EFT) QJ]. 
In EFT, only low-energy (or long-range) degrees of freedom are included explicitly, while all others are parametrized 
in terms of the most general local contact interactions. This approach relies on the fact that a low-energy probe of 
momentum k cannot resolve structures on scales smaller than 1 jk. (Note that h — c — 1 in this talk.) As a consequence, 
the influence of short-distance physics on low-energy observables can be captured in a few low-energy constants. The 
EFT describes universal low-energy physics independent of detailed assumptions about the short-distance dynamics. 

All low-energy observables can be calculated in a controlled expansion in powers of kl, where I is the characteristic 
low-energy length scale of the system. Error estimates can be obtained from the higher orders in the expansion in kl. 
The size of / depends on the system under consideration: for particles interacting via a finite range potential, e.g., I is 
given by the range of the potential. For the examples discussed in this talk, I is of the order of the effective range r e . 

We will focus on applications of EFT to few-body systems with large S-wave scattering length \a\ ^> I- For a generic 
system, the scattering length is of the same order of magnitude as the low-energy length scale I. Only a very specific 
choice of the parameters in the underlying theory (a so-called fine tuning) will generate a large scattering length a. 
Nevertheless, systems with large a can be found in many areas of physics. The fine tuning can be accidental or it can 
be controlled by an external parameter. 

Examples with an accidental fine tuning are the S-wave scattering of nucleons and of 4 He atoms. For S-wave 
nucleon-nucleon scattering, the scattering lengths and effective ranges are a s — —23.76 fm and r s = 2.75 fm in the 
spin-singlet channel and a t = 5.42 fm and r t = 1.76 fm in the spin-triplet channel. The effective ranges are both 
comparable to the characteristic low-energy length scale l K — l/m K 1.4 fm which is given by the inverse pion 
mass. However the scattering length a s is much larger than l n and a t is at least significantly larger. For 4 He atoms, 
the scattering length a w 104 A is more than a factor 10 larger than the typical low-energy length scale / ss 5 A given 
by the van der Waals interaction. The scattering length of alkali atoms close to a Feshbach resonance can be tuned 
experimentally by adjusting the external magnetic field. In Fig. ^ we show the scattering length of 85 Rb atoms as 
a function of the magnetic field near the Feshbach resonance at B = 155 G. Sufficiently close to the resonance, the 
scattering length can be made arbitrarily large. 

At very low energies all these systems behave similarly and show universal properties associated with the large 
scattering length. A simple example in the two-body sector for a > is the existence of a shallow molecule (the 
dimer) with binding energy B^ = 1/ (ma 2 ) . A particularly remarkable example in the three-body sector is the existence 
of shallow three-body bound states (Efimov states) with universal properties |2]. Efimov states can cause dramatic 
dependence of scattering observables on a and on the energy j^B]- 

We are interested in processes with typical momenta k ~ 1/M- In this case, the EFT expansion becomes an 
expansion in powers of l/\a\. In leading order, we can simply set / — and three-body observables are determined by 
a and by a three-body parameter that also determines the Efimov spectrum. The higher-order corrections in l/\a\ can 




FIGURE 1. The scattering length a of Rb atoms (in atomic units) as a function of the magnetic field close to the Feshbach 
resonance at B = 155 G (solid line). The dashed line indicates the off-resonant value of a, while the vertical dotted line gives the 
resonance position. 



be calculated but will not be discussed in detail in this talk. 

In the next section, we will briefly review the EFT for three-body systems with large scattering length a. In the 
following sections, we will discuss some applications in nuclear and atomic physics. In particular, we will discuss the 
possibility of an infrared limit cycle in QCD. Recent extensions of the EFT approach to the four-body system and 
Af-boson droplets in two spatial dimensions will also be addressed. For a more detailed review, see Ref. [4]. 



THREE-BODY SYSTEM WITH LARGE SCATTERING LENGTH 

We start with with a two-body system of nonrelativistic bosons with large S-wave scattering length a and mass m. The 
generalization to fermions is straightforward but will not be discussed in detail. In the following, we will refer to the 
bosons simply as atoms. At sufficiently low energies, the most general Lagrangian may be written as: 

* = w ' i idt + £) w ~ 8 i^ f ^ 2 | (v/tv/)3 + • ■ ■ ' (1) 

where the dots represent higher-order derivative terms which are suppressed at low energies. The strength of the 
two-body interaction g2 is determined by the scattering length a, while the three-body interaction #3 depends on a 
parameter to be introduced below. For momenta k of the order of the inverse scattering length l/\a\, the problem 
is nonperturbative in ka. The full two-atom scattering amplitude can be obtained analytically by summing the so- 
called bubble diagrams with the g2 interaction term shown in Fig. [5] The #3 term does not contribute to two-body 




FIGURE 2. The bubble diagrams with the g2 interaction contributing to the two-body scattering amplitude. 

observables. The bubble diagrams are ultraviolet divergent and have to be regulated. This can be done by introducing 
an ultraviolet cutoff A in the loop integrals. The atom-atom scattering amplitude /aa can be renormalized by allowing 
the coupling g2 to depend on A and choosing g2 (A) such that /aa is independent of the regulator. After renormalization, 
the resulting amplitude reproduces the leading order of the well-known effective range expansion for the scattering 
amplitude: /aa(&) = (— — ik)~ l , where the total energy is E — k 2 /m and k is the center-of-mass momentum. If 
a > 0, /aa nas a pole at k = i/a corresponding to a shallow dimer with binding energy B2 = 1/ (ma 2 ). Higher-order 
derivative interactions are perturbative and give the momentum-dependent terms in the effective range expansion. 



We now turn to the three-body system. At leading order in l/\a\, the atom-dimer scattering amplitude is given by 
the integral equation shown in Fig. [5] A solid line indicates the single-atom propagator and a double line indicates 
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FIGURE 3. The integral equation for the atom-dimer scattering amplitude. The single (double) line indicates the single-atom 
(two- atom) propagator. 

the full two-body propagator including all atom bubbles (cf. Fig.[2j. The integral equation contains contributions from 
both the two-body and the three-body interaction terms g2 and gj, respectively. The inhomogeneous term is given by 
the first two diagrams on the right-hand side: the one-atom exchange diagram and the contribution of the three-body 
interaction. The integral equation simply sums these diagrams to all orders. After projecting onto S-waves, we obtain 
the equation 

?( t ,r.E) - ^ M ( t . P;£ ) + l/ A ^"^ ; y^ £ ', (2, 

3fl % Jo - 1 1 a + V 3q 2 /4- — mE — IE 

for the off-shell atom-dimer scattering amplitude with the inhomogeneous term 

1 / p 2 + pk + k 2 -mE\ H(A) 
M(k,p;E) = ^{ P p2 J pk + k ^ mE )^. (3) 

The logarithmic term is the S-wave projected one-atom exchange, while the term proportional to H{A) comes from the 
three-body interaction. Note that we have introduced a cutoff regulator A in Eq. 0. This cutoff is required to insure 
that has a unique solution. It is convenient to express #3 as 

9 So 

g 3 = —$H(A), (4) 

which defines the function H(A) in Eq. Q. The S-wave atom-dimer scattering amplitude /ad is given by the solution 
& of Eq. © evaluated at the on-shell point: 

1 3k 2 1 

fAD{k) = — - = ^(k,k;E), where E=- f (5) 

Afcotoo — ik 4m ma^ 

is the total energy and k is the center-of-mass momentum. The three-body binding energies B3 are given by those 
values of E < for which the homogeneous version of Eq. Q has a nontrivial solution. 

All physical observables must be independent of the regulator A. Therefore, we can determine H(A) by varying A 
and demanding invariance of the low-energy three-body observables under this renormalization group transformation. 
This leads to the expression Ql^]: 

_ cos [50 ln(A/A„) + arctansp] 
cos [so In (A/ A*) — arctanso] ' 



where so = 1.00624... is a solution to the transcendental equation 
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FIGURE 4. The three-body coupling H as a function of the cutoff A for a fixed value of the three-body parameter A* . The solid 
line shows the analytical expression {6j, while the dots show results from the numerical solution of Eq. 



and A* is a three-body parameter introduced by dimensional transmutation. This parameter cannot be predicted by 
the EFT and must be determined from a three-body observable. The dependence of the three-body coupling H on 
the cutoff A is shown in Fig.@]for a fixed value of the three-body parameter A*. The solid line shows the analytical 
expression (|SJl, while the dots show results from the numerical solution of Eq. (0- A good agreement between both 
methods is observed, indicating that the renormalization of Eq. (0 is well under control. Note that H (A) is periodic and 
runs on a limit cycle. An important signature of an renormalization group (RG) limit cycle is discrete scale invariance. 
When A is increased by a factor of exp(7r/so) ~ 22.7, H(A) returns to its original value. Discrete scale invariance also 
arises in other contexts as varied as turbulence, sandpiles, earthquakes, and financial crashes |7j|. The possibility of 
renormalization group limit cycles was first pointed out by Wilson in 1971 01, but no explicit examples were known 
at the time. Recently, a number of field-theoretical models with limit cycles have been constructed HOi 1 1 CM. 1 1 ll 1 1 2B - 
Leclair, Roman, and Sierra have used the colorful phrase Russian doll renormalization group to describe RG limit 
cycles JTHI12I1 . The name refers to a traditional souvenir from Russia consisting of a set of hollow wooden dolls that 
can be nested inside each other. Similar to a limit cycle, the dolls show discrete scale invariance: the scaling factors 
between each doll and the next smaller one are all approximately equal. 

In summary, two parameters are required to specify a three-body system with large scattering length at leading order 
in l/\a\: they may be chosen as the scattering length a (or equivalently B2 if a > 0) and the three-body parameter A* 
. In the following sections, we will discuss some applications of this EFT in nuclear and atomic physics. 

UNIVERSAL PROPERTIES OF FEW-BODY SYSTEMS 

This EFT confirms and extends the universal predictions for the three-body system first derived by Efimov J^Sl- The 
best known example is the Efimov effect, the accumulation of infinitely-many three-body bound states at threshold 
as a — > ±°o 1 2]. The ratio of the binding energies of successively shallower states rapidly approaches a constant 
Aq = exp(2^/so) pa 515. 

However, universality also constrains three-body scattering observables. The atom-dimer scattering length, e.g., can 
be expressed in terms of a and A* as 13l 1 1 311 

ai2 = fl(1.46-2.15tan[i ln(flA*)+0.09])(l + 0{l/a)), a>0. (8) 

Note that the log-periodic dependence of H on A* is not an artefact of the renormalization and is also manifest in 
observables like 012- This dependence could, e.g., be tested experimentally with atoms close to a Feshbach resonance 
by varying a. Similar expressions can be obtained for other three-body observables such as scattering phase shifts 



as well as rates for three-body recombination and dimer relaxation (including the effects of deeply-bound two-body 
states) full. 

Since up to corrections of order l/\a\, low-energy three-body observables depend on a and A* only, they obey 
non-trivial scaling relations. If dimensionless combinations of such observables are plotted against each other, they 
must fall close to a line parametrized by A* An example of a scaling function relating 4 He trimer ground and 

excited state energies B^ and is shown in the left panel of Fig. [5] (A related scaling function was obtained in 
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FIGURE 5. The scaling functions relating the 4 He trimer ground and excited state energies (left panel) and the 4 He trimer ground 
state and the atom-dimer scattering length (right panel). The data points are calculations using various methods and 4 He potentials 
(see Ref. Hal and references therein). Note that as = 1/ i/inB^. 

Ref. Q .) The data points show calculations using various approaches and 4 He potentials. Since different potentials 
have approximately the same scattering length but include different short-distance physics, different points on this 
line correspond to different values of A*. The small deviations of the potential model calculations are mainly due to 
effective range effects. They are of the order r e /a « 10% and can be calculated at next-to-leading order in EFT. The 
calculation corresponding to the data point far off the universal curve can easily be identified as problematic since the 
deviation from universality by far exceeds the expected 10%. The right panel shows the scaling function relating the 
4 He trimer ground state energy B$ and the atom-dimer scattering length a\%. A similar scaling relation is observed 
in nuclear physics between the spin-doublet neutron-deuteron scattering length and the triton binding energy and is 
known as the Phillips line 11711 . 

Recently, we have extended the effective theory for large scattering length to the four-body system II 811 . It is 
advantageous in this case to use an effective quantum mechanics framework, since one can start directly from the well- 
known Yakubovsky equations. We have shown that no four-body parameter enters at leading order in l/\a\. Therefore 
renormalization of the three-body system automatically guarantees renormalization of the four-body system. As a 
consequence, there are universal scaling functions relating three- and four-body observables as well. As an example, 
we show the correlation between the trimer and tetramer binding energies in Fig.|S] The left panel shows the correlation 
between the ground state energies while the right panel shows the correlation between the excited state energies. The 
data points are calculations using various methods and 4 He potentials (see Ref. 1 19] and references therein). We 
conclude that universality is well satisfied in the three- and four-body systems of 4 He atoms. 

The existence of these scaling functions is a universal feature of systems with large scattering length and is 
independent of the details of the short-distance physics. Similar correlations between few-body observables appear 
for example also in nuclear systems. The correlation between the triton and a particle binding energies (the Tjon 
line [20]) and the Phillips line can be explained using a similar effective theory The spin-singlet and 

spin-triplet scattering lengths a s and a, for nucleons are both significantly larger than the range of the nuclear force. 
This observation can be used as the basis for an EFT approach to the few-nucleon problem in which the nuclear forces 
are approximated by contact interactions with strengths adjusted to reproduce the scattering lengths a s and a t l23ll . 
(For an earlier related approach, see Ref. This EFT does not contain explicit pion degrees of freedom and is 

sometimes referred to as pionless EFT. The EFT involves an isospin doublet N of Pauli fields with two independent 
two-body contact interactions: A' GjN c 7V ct OiN and TkN c N c ^ TkN, where N c — O^A* . Renormalization in the two- 
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FIGURE 6. The scaling function relating the 4 He trimer and tetramer ground state (left panel) and excited state (right panel) 
energies. The data points are calculations using various methods and 4 He potentials (see Ref. 1 19] and references therein). 



body sector requires the two coupling constants to be adjusted as functions of A to obtain the correct values of a s and 
a, . Renormalization in the three-body sector requires a three-body contact interaction A^ GiN c N cr OjNN 1 ' (7; OjN with a 
coupling constant proportional to (|6j I23II . Thus the renormalization involves an ultraviolet limit cycle. The scaling- 
violation parameter A* can be determined by using the triton binding energy B t as input. Effective range effects and 
other higher order corrections can be treated in perturbation theory II25I1 . 
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FIGURE 7. The correlation between the binding energies of the triton and the a particle (the Tjon line). The lower (upper) 
line shows our leading order result using a s and (a s and a t ) as two-body input. The grey circles and triangles show various 
calculations using phenomenological potentials 1 26]. The squares show the results of chiral EFT at NLO for different cutoffs while 
the diamond shows the N 2 LO result 1 271 12811 . The cross shows the experimental point. 

In Fig.0 we show the result for the Tjon line with a s and Bj as input (lower line) and a s and a t as two-body input 
(upper line). Both lines generate a band that gives a naive estimate of higher order corrections in £/\a\. We also show 
some calculations using phenomenological potentials 1 26] and a chiral EFT potential with explicit pions 1 271 12811 . All 
calculations with interactions that give a large scattering length must lie close to this line. Different short-distance 
physics and/or cutoff dependence should only move the results along the Tjon line. This can for example be observed 
in the NLO results with the chiral potential indicated by the squares in Fig.Qor in the few -body calculations with the 
low-momentum AW potential Viow-t carried out in Ref. 12911 . The Viow-zt potential is obtained from phenomenological 
AW interactions by intergrating out high-momentum modes above a cutoff A but leaving two-body observables (such 



as the large scattering lengths) unchanged. The results of few-body calculations with Viow-^ are not independent of A 
but lie all close to the Tjon line (cf. Fig. 2 in Ref. 12910 . 

If we take the triton binding energy B t = 8.48 MeV as input, we can use the curves in Fig.0to predict the binding 
energy of the a particle. If the spin-singlet and spin-triplet scattering lengths a s and a, are used as the two-body 
input, we find B a = 29.5 MeV. If the deuteron binding energy is used as input instead of a,, we obtain B a = 26.9 
MeV. This variation is consistent with the expected 30% accuracy of a leading order calculation in £/\a\. Our results 
agree with the (Coulomb corrected) experimental value B^ p — (29.0 ±0.1) MeV to within 10%. We conclude that 
the ground state energy of the a particle can be described by an effective theory with short-range interactions only. 
Of course, a more refined analysis should also include Coulomb corrections. However, the size of these corrections is 
expected to be smaller than the NLO contribution from the two-body effective ranges. Furthermore, corrections due to 
isospin violation will appear naturally at higher order within the effective theory. 



AN INFRARED LIMIT CYCLE IN QCD 

The low-energy few-nucleon problem can also be described by an EFT that includes explicit pion fields. Such an 
EFT has been used to extrapolate nuclear forces to the chiral limit of QCD in which the pion is exactly massless 
J30L I3H I32L l33ll . In Fig. [8] we show the results for the extrapolation of the inverse scattering lengths \/a t and l/a s 
as functions of m K from the calculation by Epelbaum et al. 13311 - The extrapolation to larger values of m K predicts 
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FIGURE 8. The inverse scattering lengths I /a, (left panel) and \/a s (right panel) as functions of as predicted by the EFT 
with pions of Ref. 1 33]. 

that a t diverges and the deuteron becomes unbound at a critical value in the range 170 MeV < m K < 210 MeV. It is 
also predicted that a s is likely to diverge and the spin-singlet deuteron become bound at some critical value of m K not 
much larger than 150 MeV. Both critical values are close to the physical value m n = 135 MeV. Beane and Savage have 
argued that the errors in the extrapolation to the chiral limit are larger than estimated in Ref. 1 33], but their errors agree 
for the small extrapolation to the critical values 

Chiral extrapolations can also be calculated using the EFT without pions The inputs required are the chiral 
extrapolations a s (m K ), a,(m K ), and B t (m K ), which can be calculated using an EFT with pions. As an illustration, 
we take the central values of the error bands for the inverse scattering lengths 1 /a s (m n ) and l/a,(m K ) from the chiral 
extrapolation in Ref. 1 33 ] . Since the chiral extrapolation of the triton binding energy B t (m K ) has not yet been calculated 
and since A* should vary smoothly with m n , we approximate it by its physical value A* = 189 MeV for m n = 138 
MeV. In Fig.|5] we show the resulting three-body spectrum in the triton channel as a function of m K . Near m n » 175 
MeV where the decay threshold comes closest to K = 0, an excited state of the triton appears. This excited state is a 
hint that the system is very close to an infrared limit cycle. In the case illustrated by Fig.|9j the value of m K at which 
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FIGURE 9. The binding momenta K = (mS 3 ) 1 1 1 of pnn bound states as a function of m n . The circles are the triton ground state 
and excited state. The crosses are the physical binding energies of the deuteron and triton. The dashed lines are the thresholds for 
decay into a nucleon plus a deuteron (left curve) or a spin-singlet deuteron (right curve). 



a t diverges is larger than that at which a s diverges. If they both diverged at the same value of m n , there would be an 
exact infrared limit cycle. 

We conjecture that QCD can be tuned to this infrared limit cycle by adjusting the up and down quark masses m„ 
and i7id 13411 . As illustrated in Fig. [9] the tuning of m n , which corresponds to m u + m c i, is likely to bring the system 
close enough to the infrared limit cycle for the triton to have one excited state. We conjecture that by adjusting the 
two parameters m u and m c i to critical values, one can make a t and a s diverge simultaneously. At this critical point, the 
deuteron and spin-singlet deuteron would both have zero binding energy and the Efimov effect |2] would occur: the 
triton would have infinitely-many increasingly-shallow excited states. The ratio of the binding energies of successively 
shallower states would rapidly approach a constant Xq close to 515. The limit cycle would be manifest in the Efimov 
effect for the triton. This argument requires that QCD with the appropriate combination of up and down quark masses 
is in the same universality class as the three-body system with large scattering length 1 35 ] . It would be interesting to 
demonstrate the existence of this infrared RG limit cycle in QCD using Lattice QCD and EFT. 



N-BOSON DROPLETS IN TWO DIMENSIONS 

In this section, we consider the universal properties of weakly interacting bosons with large scattering length (or 
equivalently a shallow dimer state) in two spatial dimensions (2D) |36]. The 2D case is very different from the three- 
dimensional case discussed above. In particular, there is no Efimov effect or Thomas collapse and no three-body 
interaction is required at leading order. We consider self -bound droplets of N(^> 1) bosons interacting weakly via an 
attractive, short-ranged pair potential. (This corresponds to Eq. Q with #3 set to zero.) Our analysis relies strongly on 
the property of asymptotic freedom of 2D bosons with an attractive interaction. 

In 2D, any attractive potential has at least one bound state. For the potential — g8 2 (r) with small g > 0, there is 
exactly one bound state with an exponentially small binding energy, 

B 2 ~A 2 exp(-47ryg), (9) 

where A is the ultraviolet momentum cutoff (which is the inverse of the range of the potential). Asymptotic freedom 
provides an elegant way to understand this result. In 2D nonrelativistic theory, the four-boson interaction g(\j/ Jf \j/) 2 is 
marginal. The coupling runs logarithmically with the length scale R, and the running can be found by performing the 
standard renormalization group (RG) procedure. For g > 0, the coupling grows in the infrared, in a manner similar to 



the QCD coupling ll37t l38ll . The dependence of the coupling on the length scale R is given by 



1 1 o i 1 ~' 

\n(A 2 R 2 ) 

g 4tt 



(10) 



-1/2 

so the coupling becomes large when R is comparable to the size of the two body bound state B 2 . This is in essence 
the phenomenon of dimensional transmutation: a dynamical scale is generated by the coupling constant and the cutoff 
scale. It is natural, then, that Z?2 is the only physical energy scale in the problem: the binding energy of three-particle, 
four-particle,... bound states are proportional to B 2 . (Remember that in 2D there is no Efimov effect and no three-body 
parameter is required.) The /V-particle binding energy Bn, however, can be very different from Z?2 if N is parametrically 
large. We use the variational method to calculate the size of the bound state. For a cluster of a large number of bosons, 
one can apply classical field theory. We thus have to minimize the expectation value of the Hamiltonian with respect to 
all field configurations \jf(r) satisfying the constraint /V = J d 2 r l/r y/ . This is equivalent to a Hartree calculation with 
the running coupling constant g(R) instead of the bare one. In the limit of a large number /V of particles in the droplet, 
some exact predictions can be obtained. 

The system possesses surprising universal properties. Namely, if one denotes the size of the /V-body droplet as R^, 
then at large /V and in the limit of zero range of the interaction potential 1 3(1 : 

-jr 1 ^ 0.3417, -^ ±i ~ 8.567, /V>1. (11) 
Rn Bjm 

The size of the bound state decreases exponentially with /V: adding a boson into an existing /V-boson droplet reduces 
the size of the droplet by almost a factor of three. Correspondingly, the binding energy of /V bosons Bn increases 
exponentially with N. This implies that the energy required to remove one particle from a /V-body bound state (the 
analog of the nucleon separation energy for nuclei) is about 88% of the total binding energy. This is in contrast to 
most other physical systems, where separating one particle costs much less energy than the total binding energy, 
provided the number of particles in the bound state is large. Equations il Q are accurate to leading order in l/N but 
the 1 //V-corrections are calculable. 

For the universal predictions il Tj to apply in realistic systems with finite-range interactions, the /V-body bound 
states need to be sufficiently shallow and hence have a size Rn large compared to the natural low-energy length scale 
/. Depending on the physical system, / can be the van der Waals length, the range of the potential or some other scale. 
As a consequence, Eqs. (lilt are valid in realistic systems for /V large, but below a critical value, 

1 «/V«/V crit « 0.931 \n{R 2 /l) + 0{\). (12) 

AtN = /Vcrit the size of the droplet is comparable to I and universality is lost. If there is an exponentially large separation 
between R2 and /, then N cr it is much larger than one and the condition H2\ can be satisfied. 

We can compare our prediction with exact few -body calculations for /V = 3,4. While the 1/iV-corrections are 
expected to be relatively large in this case, we can estimate how the universal result for B^ + \/Bn is approached. 

The three-body system for a zero-range potential in 2D has exactly two bound states: the ground state with of' = 

16.522688(1)B 2 and one excited state with = 1. 2704091 (l)B 2 JHEIS. Similarly, the four-body system for 

a zero-range potential in 2D has two bound states: the ground state with B^ = 197.3(1)Z?2 and one excited state 

with B^ =25.5(1)52 lEUl . The prediction il It applies to the ground state energies B^ 1 and . The exact few- 
body results for /V = 3,4 are compared to the asymptotic prediction for Bn/Bn-i indicated by the dot-dashed line in 

Fig. [TO] The ratio B^ /B2 ~ 16.5 is almost twice as large as the asymptotic value (lilt , while the ratio B3 /B4 « 1 1 .9 
is already considerably closer. These deviations are expected for such small values of N. Note, however, that the ratio of 
the root mean square radii of the two- and three-body wave functions is 0.306 1 40], close to the asymptotic value il It . 
The dashed line gives an estimate of how the large-/V value should be approached 14 ill . This estimate assumes the 
expansion: 

b n= (cq + ^tt + ^ l + --) 8.567^, (13) 
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FIGURE 10. The ratio Bn/Bn-\ as a function of N. The dot-dashed line shows the asymptotic value of 8.567. The dashed line 
is an estimate of how this value is approached as N increases. 



The dashed line in Fig.^|was obtained by fitting the coefficients of the 1 /N 2 and 1 /N^ terms in Eq. dl4> to the data 
points for Af = 3,4. 

It would be interesting to test the universal predictions Jl Q both theoretically and experimentally for N > 4-. On the 
theoretical side, Monte Carlo techniques appear to be a promising avenue. Furthermore, the experimental realizablity 
of self -bound 2D boson systems with weak interactions should be investigated. According to the estimate in Fig. 1101 
the 1 /N corrections to Eqs. (II It are small for N > 6. Using (I12> . this requires R%fi 3> 600. We are not aware of any 
physical system that satisfies this constraint. However, such a system could possibly be realized close to a Feshbach 
resonance where R2 can be made arbitrarily large. 



SUMMARY 

We have discussed the EFT for few-body systems with short-range interactions and large scattering length a. The 
renormalization of the three-body system with large a in three spatial dimensions requires a one-parameter three-body 
force governed by a limit cycle already at leading order in the expansion in l/\a\. As a consequence, two parameters 
are required to specify a three-body system: the scattering length a (or the dimer binding energy B2) and the three- 
body parameter A*. Once these two parameters are given, the properties of the three- and four-body systems are fully 
determined at leading order in l/\a\. 

The large scattering length leads to universal properties independent of the short-distance dynamics. In particular, 
we have discussed universal expressions for three-body observables and universal scaling functions relating various 
few-body observables for both atomic and nuclear systems. A more detailed account including the effects of deeply- 
bound two-body states can be found in Ref. [4]. 

The success of the pionless EFT demonstrates that physical QCD is close to an infrared limit cycle. We have 
conjectured, that QCD could be tuned to the critical trajectory for the limit cycle by adjusting the up and down 
quark masses. The limit cycle would them be manifest in the Efimov effect for the triton |34]. It may be possible to 
demonstrate the existence of this infrared RG limit cycle in QCD using Lattice QCD and EFT. 

In two spatial dimensions, the three-body parameter A* does not enter at leading order in the expansion in l/\a\ 
and A^-body binding energies only depend on B^- The asymptotic freedom of non-relativistic bosons with attractive 
interactions in 2D leads to remarkable universal properties of A^-body droplets, such as the exponential behavior of 
binding energies and droplet sizes in Eq. (II It |36[. 

The three-body effects discussed here will also become relevant in Fermi systems with three or more spin states. 
Future challenges include universality in the A^-body problem for N > 4, effective range corrections, and a large 



scattering length in higher p artial waves. A large P-wave scattering length, for example, appears in nuclear Halo 
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